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Abstract 

Loop Quantum Cosmology (LQC), mainly due to Bojowald, is not the cosmological sector 
of Loop Quantum Gravity (LQG). Rather, LQC consists of a truncation of the phase space of 
classical General Relativity to spatially homogeneous situations which is then quantized by the 
methods of LQG. Thus, LQC is a quantum mechanical toy model (finite number of degrees of 
freedom) for LQG (a genuine QFT with an infinite number of degrees of freedom) which provides 
important consistency checks. However, it is a non trivial question whether the predictions of 
LQC are robust after switching on the inhomogeneous fluctuations present in full LQG. 

Two of the most spectacular findings of LQC are that 1. the inverse scale factor is bounded 
from above on zero volume eigenstates which hints at the avoidance of the local curvature 
singularity and 2. that the Quantum Einstein Equations are non - singular which hints at 
the avoidance of the global initial singularity. This rests on 1. a key technique developed for 
LQG and 2. the fact that there are no inhomogeneous excitations. We display the result of a 
calculation for LQG which proves that the (analogon of the) inverse scale factor, while densely 
defined, is not bounded from above on zero volume eigenstates. Thus, in full LQG, if curvature 
singularity avoidance is realized, then not in this simple way. 

In fact, it turns out that the boundedness of the inverse scale factor is neither necessary nor 
sufficient for curvature singularity avoidance and that non - singular evolution equations are 
neither necessary nor sufficient for initial singularity avoidance because none of these criteria 
are formulated in terms of observable quantities. After outlining what would be required, we 
present the results of a calculation for LQG which could be a first indication that our criteria 
at least for curvature singularity avoidance are satisfied in LQG. 

*jbrunnem@aei.mpg.de,jbrunnemann@perinieterinstitute.ca 
^thieinann@aei.mpg.de,tthiemann@perimeterinstitute.ca 
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1 Introduction 



Loop Quantum Gravity (LQG) is a candidate for a Quantum Field Theory (QFT) in four dimensions 
which achieves to unify the principles of Quantum Theory (QT) and General Relativity (GR). This 
means that LQG implements the fundamental feature of GR, its background independence, in a 
quantum setting, see e.g. ^ for books and for reviews. 

To see how radical that ambition is, recall that ordinary QFT relies on a background spacetime 
(M, qq). Here M is a differential manifold and go is a prescribed background metric thereon. Ordinary 
QFT now axiomatically assumes causality with respect to qq: Given two smeared field operators 
such that the supports of the smearing functions /, /' are spacelike separated (the non - 
spacehke geodesies with respect to go starting from supp(/) never hit supp(/') and vice versa) then 
the field operators must (anti)commute. Thus we see that the background metric go pivotally finds 
its way into the very definition of the field algebra. Without go we therefore do not even know what 
a quantum field is (see however [S] for an interesting reformulation of ordinary QFT which may get 
us rid of go) ■ 

The background go, however, is incompatible with GR, in particular Einstein's equations, which 
state that the metric g is a dynamical field which cannot be prescribed but must be determined 
dynamically according to the matter energy momentum density. Hence, ordinary QFT neglects 
backreaction effects from matter on geometry (which is appropriate in scattering situations of negible 
particle numbers). Worse, if g becomes itself a quantum field, in particular if there is no background 
go such that the fluctuations g — go are small in an appropriate sense, then ordinary QFT fails to 
be a valid description of the physical processes^ . It is precisely in those situations that a full fledged 
quantum theory of both matter and geometry must take over and it is almost granted that these 
situations arise for instance in the early universe. 

Thus, the task of a quantum gravity theory must be to define a quantum field just on a differential 
manifold M and not on a background spacetime (M, go) . This means that one has to throw many 
aspects of conventional QFT over board. LQG is trying to precisely achieve this goal. As one 
might expect, given that already classical GR is a rather difficult classical field theory, this is a very 
hard enterprise, both conceptually and technically. Despite these difficulties, LQG has made steady 
progress over the past fifteen years and we are getting closer and closer to being able to answer the 
question whether LQG really is a QFT of GR in the sense that 1. classical GR is the classical limit 
of LQG and 2. when the fluctuations of the quantum metric are small, then the framework of QFT 
on (M, go) is recovered. In other words, if these questions could be answered affirmatively, then we 
would have a viable quantum gravity theory in front of us and we could start drawing conclusions 
and physical predictions from it and hopefully compare those with experiments. 

However, even before we (hopefully) complete the construction of full LQG, it is possible and 
mandatory to test the theory by models which are conceptually and technically easier to handle and 
which still capture enough of the features of the full theory. To give an example, the hydrogen atom 
can be analyzed just using quantum mechanics for the electrons orbiting a classical nucleus and one 
gets fantastically close to the experiment. Hence, this gives us an important first idea. However, the 
full problem of the hydrogen atom must certainly take into account the full QFT of the standard 
modeP and only then one can determine the Feinstruktur of the hydrogen atom energy spectrum 

"'^The alternative, to quantize the fluctuations themselves and to construct a graviton QFT on go fails due to 1. 
the non - renormalizability of GR (see however 0] for ideas concerning a non - trivial UV fix point) and 2. in violent 
situations there may simply be no go at all such that the fluctuations are small. 

^Even today that has not yet been done fully, usually one only considers QED effects but not e.g. the QCD physics 
of the proton. 
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such as the Zitterbewegung of the electron due to the interaction with the electromagnetic field. 

It is remarkable and important that the predictions of the simple quantum mechanical compu- 
tation are confirmed by the full theory because a priori it could be that the full QED calculation 
completely changes the picture. Of course, in the case of the hydrogen atom one has access to ex- 
perimental data and so this was not expected if QED has something to do with nature. However, 
in case of absence of data there is no reason why one should trust a toy model calculation at all. 
In particular, it is unclear why a quantum mechanical calculation describing only a finite number of 
degrees of freedom will reproduce the QFT aspects of the full problem involving an infinite number 
of degrees of freedom. One hopes of course that the full problem will only yield corrections to the 
model, however, that is far from granted and must be proved rigorously. 

As far as LQG is concerned, in recent years a new class of toy models have been constructed which 
were coined Loop Quantum Cosmology (LQC). These developments are mainly due to Bojowald, see 
e.g. the beautiful review . These models consider the truncation of the infinite dimensional phase 
space of GR to spatially homogeneous situations which is finite dimensional. Such truncations have 
been considered before and are called minisuperspace models jH]. Traditionally they were quantized 
using the Schrodinger representation of the canonical commutation relations. What is new in Bo- 
jowald's work is that he employs a different representation of the canonical commutation relations. 
This might come as a surprise because the Stone von Neumann theorem gives us a uniqueness result 
concerning the representation theory of the Weyl algebra for finite dimensional phase spaces. The 
catch is that the Stone von Neumann theorem assumes, among other things, (weak) continuity of the 
Weyl operators. Bojowald drops this assumption because in full LQG the unique [7j representation 
[H] is also discontinuous, hence this representation is a much better model for LQG than the usual 
Schrodinger representation. This kind of representations were previously discussed by e.g. Thirring 
[5] and while exotic are used e.g. in solid state physics (TU]. Moreover, using those one can circumvent 
negative norm states in Maxwell theory ^1] or string theory which then does not lead to restrictions 
on the spacetime dimension ^2]- 

A second important ingredient in Bojowald's work is a key technique discovered for full LQG 
in which allows to define certain operators densely on the LQG Hilbert space although their 
classical counterpart is ill defined when the metric becomes degenerate. Prime examples for such 
classically singular situations are cosmological big bang singularities. This key technique was used in 
[T!^ in order to give, for the first time, a mathematical meaning to the Quantum Einstein Equations 
in full LQG, mathematically speaking, the so - called Wheeler - DeWitt equations could be defined 
densely on the LQG Hilbert space. Moreover, these equations can be solved explicitly by following a 
constructive algorithm J3] . This works for arbitrary matter coupling ^S] and was tested successfully 
in 2+1 gravity ^H]- For a modern version of this framework see ^7]. A simpler setting in which this 
key technique was used is the construction of the length operator in full LQG jTHj. This operator 
is also well (i.e. densely) defined, it is an unbounded but positive self - adjoint operator whose 
spectrum is entirely discrete. 

Of course, due to the complexity of the full theory, LQG cannot be solved completely. It is 
precisely the virtue of toy models such as LQC that their mathematics is comparatively simple 
so that one can focus immediately on the conceptual analysis. Bojowald therefore could go very 
far in the quantization of the cosmological models although some important aspects such as the 
physical inner product are still missing. The most spectacular finding is that in LQC the operator 
corresponding to the inverse scale factor(s)^ is (are) bounded from above [^|2n|- This hints at the 
following: Recall that in classical GR one characterizes singularities by global and local criteria |21j . 

■^One for the isotropic models, three for the diagonal models. 
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The global criterion is the causal geodesic incompleteness of a (inextendible) spacetime. Equivalently, 
the Einstein equations on globally hyperbolic spactimes with manifolds diffeomorphic to R x a cannot 
be solved for all values of t G 11. The local criterion is defined as the divergence of polynomial scalars 
built out of the Riemann tensor and is used to characterize a global singularity in more detail. The 
local criterion by itself is unsatisfactory because a spacetime can have a global singularity while all 
curvature scalars are regular. Notice that both criteria are to be measured by an obeserver and in 
this sense are to be formulated in terms of observables (measurable quantities). 
In classical cosmology both types of singularities are present. Since in cosmology curvature scalars 
are given by powers of the inverse scale factor, its boundedness hints at the absence of the local 
singularity in LQC in the sense that with respect to any kinematical state the curvature expectation 
value can never diverge. Moreover, the truncated Quantum Einstein Equations (Wheeler DeWitt 
Equation) do not become singular at the classical Big Bang and one therefore can go before the Big 
Bang. This hints at global singularity avoidance in LQC. Further consequences of the boundedness of 
the inverse scale factor are that the quantum geometry behaves very non - classically at the classical 
singularity implying effective modifications of the Friedmann equations which then drive geometric 
inflation without the necessity of fine tuning the inflaton potential |22] ■ These results related to 
the boundedness of the inverse scale factor in LQC are so promising and spectacular that it was 
extensively highlighted in the physics literature and public press. In particular, the boundedness was 
mentioned as crucial to establish the non singularity of the quantum evolution equations, see e.g. 

In it was shown that it is precisely those two ingredients imported from full LQG into 
LQC, 1. the non - standard representation of the canonical commutation relations and 2. the key 
technique [T^ to densely define classically degenerate operators, which are responsible for the absence 
of the local singularity in LQC and similar models inspired more by a metric than connection based 
approach. However, as outlined above, on general grounds one must critically examine the stability 
of such model calculations when one reintroduces the field theory degrees of freedom. This has 
been stressed before for geometrodynamics in [22] where one generically finds that the model even 
quantitatively has large deviations from the full theory or more complete theory. More recently also 
LQC has been examined critically [211 concerning the issue of the semiclassical limit but here the 
results seem to be inconclusive so far. 

The main purpose of this paper is to contribute to the stability analysis of LQC versus LQG. 

The paper is organized as follows: 

In section two we summarize and compare the main features of LQG and LQC for the unfamil- 
iar reader. Experts can safely skip this section. 

In section three we report the results of a computation carried out in our companion paper [2H] 
within full LQG which is analogous to the one carried out in LQC. In LQC this calculation proves 
the boundedness of the inverse scale factor. More precisely, we ask, as in LQC, whether the part 
of the energy momentum operator which probes the quantum metric is bounded from above on the 
kinematical Hilbert space of full LQG, at least at the classical singularity, that is, when the volume 
of space vanishes. This precisely mirrors the calculation of [12]. The result is negative: In full 
LQG the (analogon of) the inverse scale factor is unbounded from above even when the 
quantum volume vanishes. This is true even when we compute its norm with respect to a state of 
zero volume which is homogeneous and isotropic on large scales. This proves that the boundedness 
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of the inverse scale factor in isotropic and homogeneous LQC does not extend to the full theory even 
when restricting LQG to those states which one would use to describe a maximally homogeneous 
and isotropic situation (modulo fluctuations). 

It is true that in non - isotropic, homogeneous models the inverse scale factor is also unbounded 
on the full Hilbert space. However, in those more general models it is still bounded on zero volume 
eigenstates which again seems to indicate that the local singularity is evaded in LQC. Our result 
is stronger: No matter whether homogeneous and/or isotropic, spherically symmetric, cylindrically 
symmetric etc., the inverse scale factor is unbounded on zero volume eigenstates in full LQG and 
on the full Hilbert space anyway. The reason for this is that LQG, in contrast to LQC, admits 
inhomogeneous, microscopical excitations which precisely account for the unboundedness. We will 
see this explicitly when we construct the zero volume eigenstates on which the inverse scale factor 
has arbitrarily large norm. Notice that meanwhile also inhomogeneous models such as spherically 
symmetric or cylindrically symmetric ones have been quantized by LQG methods in which the inverse 
scale factor apparently also stays bounded [2Zj. While these do have local, inhomogeneous degrees of 
freedom (in the spherically symmetric case without matter, at least before solving the constraints), 
as our calculation reveals, the inhomogeneous excitations of full LQG are more general than in those 
models and thus give unboundedness even on states of LQG describing a sector appropriate for those 
more general models. 

However, this does not mean that LQG does not predict the absence of the local initial singu- 
larity. Namely, in section four we remark that an operator may be unbounded but still it may be 
bounded when restricted to a subspace of the Hilbert space. Thus, the boundedness of the inverse 
scale factor is not a necessary criterion for local singularity avoidance. In our case we are interested in 
a sector which describes a collapsing universe which is homogeneous and isotropic at large scales. We 
perform a corresponding calculation which may be taken as an indication that the expectation value 
of the analogon of the inverse scale factor in LQG, with respect to kinematical^, coherent states [33] 
peaked on homogeneous and isotropic initail data, is bounded from above at the Big Bang^. 
We did this for scalar matter only but qualitatively nothing changes in more general situations. Also 
our result is completely general: We derive a bound of the inverse scale factor on coherent states 
peaked on an arbitrary point in phase space. This formula can then be specialized also to inhomo- 
geneous models and for the Bianchi I (Kasner) case it is easy to see that one gets boundedness as well. 

In section five we stress that the results of section four are promising but inconclusive: Unfortunately, 
the boundedness of the inverse scale factor is also not a sufficient criterion for local singularity avoid- 
ance. This has to do with the fact that both singularity criteria fundamentally have to be discussed 
at the level of physical observables and physical states. They must be be disussed separately because 
the quantum theory could avoid one of them but not the other. 

Concerning the local singularity, the inverse scale factor is not gauge invariant and therefore not an 
observable. It can be turned into an observable using the technique of partial observables (23 EO] but 
then its spectrum on the physical Hilbert space can differ drastically from the kinematical spectrum 
pij . Hence, reliable statements about the local singularity can only be obtained on the physical 
Hilbert space. 

As far as the global singularity is concerned, physical states automatically solve the quantum Ein- 

^Kinematical states, in contrast to physical states, do not solve the Quantum Einstein Equations. Non - gauge 
invariant operators such as the inverse scale factor can only be probed on the kinematical Hilbert space. 

^Notice that coherent states peaked on homogeneous and isotropic initial data still have inhomogeneous and 
anisotropic excitations, however, they are small. 
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stein equations so there seems to be no sign of the global singularity at the level of physical states, 
by construction^. More precisely, in order to solve the Quantum Einstein Equations, one takes an 
arbitrary superposition of kinematical states and then obtains a recursion relation for the coefficients. 
These have been derived explicitly for full LQG in ^1] and those of LQC ^Hl are similar but simpler. 
One of the labels of these coefficients can be identified with an unphysical time parameter. If one 
chooses that label in such a way that its range is (a discrete subset of) the full real axis, a solution 
to the Quantum Einstein Equations, that is, a physical state, describes an entire quantum universe 
(a history) and involves the full range of the unphysical time parameter. Thus there is no sign of 
the global singularity within a given solution even in the full theory. However, following ^H] one may 
argue that the existence of a global singularity could manifest itself in the fact that the recursion 
relations break down at zero volume. Such a breakdown can manifest itself in the following two ways: 
A. 

It can lead to a restriction or additional consistency condition on the freedom in the choice of the 
initial data (coefficients) that parameterize the solution space. This makes the solution space smaller 
than it would be otherwise. In particular, if the number of additional conditions is larger than the 
free initial parameters of the solution, the solution space would be empty. 
B. 

It can happen that certain coefficients remain undetermined by the recursion relation and lead to 
an indeterministic quantum evolution (given initial coefficients, we cannot construct the full solution 
without making additional choices). This makes the solution space larger than it would be otherwise. 
Notice that this undeterministic quantum evolution corresponds to an undeterministic classical evo- 
lution (initial value problem). However, by itself such ambiguities are not physically relevant because 
this is not an evolution of physical quantities. Even classically it is an evolution with respect to an 
unphysical time of unphysical quantities. In terms of physical (gauge invariant) quantities there is 
no unphysical evolution at all by definition because Dirac observables Poisson commute with the 
Hamiltonian constraint. Hence, what we actually have to do in order to observe a possible break- 
down of evolution is to consider an evolution with respect to a physical Hamiltonian and of physical 
observables. We will discuss this at length in section five. 

We now point out that absence of either effect A or B is neither a sufficient nor a necessary criterion 
for global singularity avoidance. 

To see that it is not necessary, notice that the presence of either effect indeed affects the size of the 
set of formal solutions to the evolution equations. However, even if such effects exist, it may still 
be the case that the physical Hilbert space is large enough in order to accomodate all semiclassical 
physical states which describe all classical spacetimes (at large scales). To see that it is not sufficient, 
notice that not all formal solutions to the Quantum Einstein Equations correspond to physical states 
because they might not be normalizable with respect to the physical inner product or they might 
have zero norm. In both cases we must drop them from the set of physical states. Hence, even 
if effects A or B are absent, it could turn out, in the worst case, that most of the solutions that 
one finds do not lie in the physical Hilbert space and then again one would be confronted with a 
global singularity because there could be not even be a semiclassical sector describing cosmology. 
Furthermore, as argued above, global singularity avoidance must be formulated in terms of physical 
observables which has not been done yet. Therefore, the breakdown of the recursion or the absence 
thereof is inconclusive for the presence or absence of global singularities unless one knows the physical 
Hilbert space and the corresponding Dirac observables. 

We conclude: In order to decide on the presence or absence of both the local and global 

^Of course, the inverse scale factor must be densely defined to even define the Quantum Einstein Equations but 
that has been established for LQG in jl3[ll5j . 
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singularity, detailed knowledge of the physical Hibert space is a necessary prerequisite.. 

This Hilbert space is known to exist both in LQG and LQC, following for instance ^ITj. However, the 
corresponding physical inner product is rather difficult to construct explicitly, for any of the current 
versions of the Hamiltonian constraint, which is why a definitive conclusion on the singularity issue 
is not available at the moment. 

To improve on this, we first of all broadly outline which steps should be fundamentally performed 
in order to prove the absence of both the global and local Big Bang singularity in full LQG. We give 
both a precise and an approximate scheme. The ideal (precise) procedure involves the exact knowl- 
edge of the Hilbert space of physical states and a sufficiently large set of Dirac observables (gauge 
invariant quantities) which, as we said, are not yet available neither in full LQG nor in LQC in a form 
explicit enough. The approximate scheme stays within the kinematical Hilbert space over which one 
has excellent control. The calculation performed in section four can be viewed as a simplified version 
of the approximate scheme. The application of either scheme is left for future research and again 
LQC would be an ideal testing ground for this scheme. 

We conclude in section six summarizing the status of the absence of the initial singularity in full 
LQG. Notice that qualitatively all we have said is also applicable to other singularities such as those 
connected with black holes In fact, most considerations within LQG concerning black holes 

are based on minisuperspace calculations which exploit the fact that the interior of a black hole can 
be mapped into a homogeneous Kantowski - Sachs model so that one can transfer, almost literally, 
the LQC results. However, our remarks prevail: No reliable conclusions can be drawn before one 
works with physical operators and physical states of the full theory. This has also been stressed in j3S]- 

Notice that this paper is not meant as a criticism of LQC by itself, in the contrary, LQC has 
provided us with many fruitful new conceptual ideas and one would hope that its results hold also 
LQG, albeit in a technically different realization. However, we want to draw attention to the fact 
that LQC is far from being full LQG and that therefore caution should be applied when trying to 
extrapolate LQC results to LQG. This extends also to more general reduced models such as the 
spherically symmetric or cylindrically symmetric ones. 

2 Relation between LQG and LQC 

For the benefit of the reader unfamiliar with the basic structure of either LQG or LQC we here 
summarize elements of both frameworks and compare them. The presentation will be oversimplified 
and the mathematical details will be discared. As we will see, LQC does not describe a stable 
subsector of LQG, rather it is a truncation in which all but finitely many degrees of freedom are 
cancelled by hand. Experts can safely skip this section. 

2.1 Elements of LQG 

The classical starting point of LQG is a Hamiltonian formulation of the Einstein Hilbert action, 
however, written in unusual variables. Instead of using the three metric qab on the spatial slices a of 
a foliation of the four manifold M = 'Kx a and the extrinsic curvature Kab as the canonical variables 
as advocated by Arnowitt, Deser and Misner, one uses variables that are more famililar from the 
canonical formulation of Yang Mills theories. They consist of an SU{2) connection and an electric 
field E"-. Here a,b,c, .. = 1, 2, 3 are tenorial indices and j, k,l, .. = 1, 2, 3 are su{2) indices. The origin 



7 



of SU{2) is that SU{2) is the universal covering group of SO (3) (necessary for couphng of spinorial 
matter) and 5*0(3) appears naturally when we write the three metric in terms of cotriads (or frame 
fields) e{, that is, gab = ^i^b^jk- Indeed, the cotriads are determined by Qab only up to an 5*0(3) 
rotation and this is how SU (2) finds its way into this formulation as a gauge group. 
The relation between {A, E) and (g, K) is 

Ai = Ti + pKabe'^, E^ = ./d^e'] (2.1) 

Here F is the spin connection (a certain function of and its first spatial derivatives), the triad e'j is 
the inverse of e{, i.e. e'^el = 6^, e^e^ = Sj and /3 > is a parameter called the Immirzi parameter. 
That {A, E) form a canonical pair means that 

{E'^{x),At{y)} = SvrG/? S^, S{x,y) (2.2) 

all others vanishing. The phase space coordinatized by {A, E) is subject to the usual spatial diffeo- 
morphism constraint which in these variables can be written as 

O, = Fi, E] where F^, = d^A^, - d^Ai + e.^iA^X (2.3) 

is the curvature of A, the Hamiltonian constraint 

^ _ ^MFL - (1 + ^')ejmnKTKnE^,E^ 

where = A{ — r^(£') and the additional Gauss constraint 

Cj = daE;; + ejkiA'.E^ (2.5) 

which gets us rid of the SU{2) degrees of freedom. We see that GR can be cast into the form of 
a SU{2) Yang Mills theory (which is also subject to Cj = 0) with one important difference: The 
background dependent Yang Mills Hamiltonian density 

H = -j^= [E^El + Bpl]P'^ (2.6) 
v/det(gO) 

where = e°'^''FlJ2 is the magnetic field and g° is the spatial projection of the background metric 
5f°, is replaced by the background independent constraints O^, C. 

If there is matter present such as a scalar field then the constraints Oq, C are augmented by 
additional terms 

Or"^^ = 7r0,, and C-^'^^ = ^^^^^2^=^ + V\d^\V{<P) (2.7) 

2 ^|det(E)| 

where vr is the momentum conjugate to (j). 

Canonical quantization now, very roughly, means to find a Hilbert space Ti-Kin on which the 
canonical Poisson brackets ()2.2|) are implemented as canonical commutation relations and to impose 
the quantum constraints, e.g. Ctp = 0. The resulting space of solutions is the physical Hilbert space 
on which the (Dirac) observables, that is, the gauge invariant functions F on phase space satisfying 
{Ca{x), F} = {C{x), F} = {Cj{x), F} = act as self adjoint operators F. This simple recipe has to 
be refined in QFT because of the singular short distance behaviour of the fields A, E. In order to 
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deal with this problem one has to smear the fields. There are many ways to smear the fields, each 
of which gives rise to a different algebra 21 of elementary observables''. Guided by experience with 
canonical lattice gauge field theory in LQG one considers the following objects: 

1. For each curve e in cr the holonomy of the connection A along e, that is, 

r ^ pi 1*1 nl 

A{e) := V exp( / A) := I2 + / dti / dt2 .. / dt^ (2.8) 

J e Jo Jti J t„_i 

where e : [0,1] —>■ e;t e{t) is a parameterization of the curve and A{t) := Ai^{e{t))e°'{t)Tj with 
Tj = —iaj/2 a basis of su{2) {aj are the Pauli matrices). 

2. For each surface S and each su{2) valued function x f^{x) the electric fiux of E through S 

EjiS) = I d'u r{Xs{u)) eatcE^iXsiu)) X^^„.(n) X^^A^) (2.9) 

J s 

where Xs : s C S; u = {u^,u'^) ^— Xs{u) is a parameterization of S. 

Using (j2.2p one can show that the Poisson brackets* of the holonomies and fiuxes close among each 
other so that they define a well defined algebra 2t. 

The kinematical Hilbert space of LQG can now be described very easily as follows: 
A union of curves e forms a graph 7. By subdividing the curves we may assume that the curves are 
mutually disjoint except for the endpoints. These curves are then called edges of 7 which intersect 
in the vertices. We will denote by -^(7), ^(7) the set of edges and vertices of a graph respectively. 
The wave functions are then of the form 

^{A)=^,{A{e,),..,A{eN)) (2.10) 

where 7 can be any graph, is the number of edges of the graph 7 and : SU{2)^ —>■ (C is a. 
complex valued function on N copies of SU{2). 
The holonomy operator acts by multiplication 

[A{^) ^]{A) := A{e) (2.11) 

while the fiuxes act by derivation 

[eJ{S) ij]{A) := th{Ef{S),m)} (2.12) 

The scalar product between functions ip, ip' defined via (j2.im) over graphs 7, 7' respectively is defined 
as follows: Take the union graph 7" := 7 U 7' and write the edges e, e' of 7,7' respectively as 
compositions of edges e" of 7", say e = e" o .. o e". Now use the fact that the holonomy factorizes, 
e.g. A{e) = A{e'l)..A{e'^) to write ipjip' as functions over 7". Then 



<ilj,ij'>:= dfinih) .. dj^nifiN) i^'y"ihi,..,hN) 'ilj'^„{hi,..,hN) (2.13) 

where N is the number of edges of 7" and fin is the Haar measure on SU (2). Specifically, if we write 
h = cos(x)l2 + Tjn^{9, ifi) sin(x) with ^ ^ [0, tt) and (p G [0, 27r) where 

^They are called elementary because one can write every function on phase space in terms of (limits of) them. 
They are not physically observable because they are not gauge invariant. 
®More precisely the corresponding Hamiltonian vector fields. 
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n = (sin(6') cos((y9), sin(6') sin(y9), cos(6')) is the usual unit vector on the sphere then duuih) = 
csin^(x) sm{9)dx dO dip and c is a normahzation constant fixed by requiring that iih{.SU{2)) = 1. 

It turns out that the Hilbert space Ti-xin has a convenient orthonormal basis consisting of the 
spin network functions (SNF) T j f. For the purpose of this paper it is enough to say that they are 
specific functions of the type ()2.1()j) which are polynomials in the holonomies A(e), e E E('j). They 
are labelled by a graph 7, a collection j := {je}ee-B(7) of spin quantum numbers je, 2je = 1, 2, 3, .., 
one for each edge, and a collection of other quantum numbers / = {Iv}vi^v('y) called intertwiners, one 
for each vertex. Since the label 7 is continuous, the basis is uncountable and the Hilbert space is 
therefore not separable. We will use the short hand s = (7(5), j(s), /(s)) for a spin network label. 

Let us now turn to the quantum constraints. The Gauss constraint imposes certain restrictions on 
the quantum numbers / of the SNF and is easily taken care of. The spatial diffeomorphism constraint 
enforces that we should identify all spin network labels s, s' if there exists a spatial diffeomorphism 
ip of a such that ^9(7(5)) = 7(5') and so is also easily taken care of. The real challenge for LQG 
(and any other quantum gravity theory) is the Hamiltonian constraint. Looking at ()2.4|) the mere 
definition of the operator C looks like a hopeless task due to the non - linearity of the constraint. It 
is here where the key technique of jTSj comes into play: 

It is easy to see that the nonlinear expression involving the electric field E can be written in the 
compact Poisson bracket form 

< = ^ ^> ^^^^^ ^ = / d'x./\d^\ (2.14) 

A/|det(£/)| J a 

is the volume^ of cr. The advantage of this way of writing ()2.14jl is that there exists a volume 
operarator corresponding to V which can be diagonalized in terms of spin network states. Thus, 
(r2.14|l can be quantized by substituting Poisson brackets by commutators divided by ih. Next, 
consider the smeared quantity C{N) = J d^xN{x)C{x) where is a test function. We may write 
this as C{N) = Cp{N)+Ck{N) where Cp, Ck respectively correspond to the terms in ()2.4p involving 
the term F^^ or K'^KJ^ respectively. Now it is easy to see that 

Ki{x)^{Ai{x),{CF{l),V}} (2.15) 

where Cf{1) := Cf{N)n=i. It follows that we know how to quantize C{N) once we know how to 
quantize CpiNy^. In order to do this it remains to write Cf{N) in terms of holonomies. This is 
achieved by recalling that for a curve e"„(t) := + tu"' and loop ax,u,v(y^) '■= ex,„(e) o ex+eu,v{^) o 
ex+e(u+v),u{^)~^ ° Gx+ev,vi^)~^ rcspcctively we have the expansions A{ex^u{^)) = I2 + eu°' Ai^{x)Tj / 2 + 
0(e2) and A(a^,„,„(e)) = I2 + e^u"'v^F^f^{x)Tj/2 + O(e^). Therefore, if we write the integral Cf{N) 
as the limit of a Riemann sum over boxes of coordinate volume then in the limit it is appropriate 
to replace by a holonomy along a loop on the boundary of the box while the A{ in the Poisson 
bracket {^^, V} may be replaced by a holonomy alon an edge on the boundary of that box. 

The details of taking that limit can be inferred from ^] • Roughly speaking, one defines (a 
spatially diffeomorphism invariant version of) the constraint on the space of solutions to the diffeo- 
morphism constraint on which the limit becomes trivial because what matters to the constraint are 
only the terms of the Riemann sum next neighbouring the vertices and those terms are diffeomorphic 
to each other when we refine the Riemann sum. 



^In the case oi k = this diverges classically. However, we only need its Poisson bracket which is well defined. 
^°The condition C{x) = for all x is equivalent to requiring C{N) = for all N because the support of N can be 
localized arbitrarily. 
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2.2 Elements of LQC and Comparison with LQG 



LQC, by definition, is the canonical quantization of tlie cosmological truncation of classical GR by 
the methods of LQG^^. Hence one first specializes the functions A{{x), Ej{x) to the spatially homo- 
geneous case, plugs those into the formulae for the holonomies and fluxes and then computes the 
corresponding algebra ^cosmo- That algebra is of course tremendously smaller than the full algebra 21 
but one quatizes it in analogy to LQG, that is, one chooses a Hilbert space representation which is as 
similar as possible to the one of LQG. Finally, one plugs the simplified expression for {A, E) into the 
constraints and quantizes them as in LQG, in particular one uses the technique fl2.14|) which enables 
to avoid denominators (which are singular at the classical singularity) at the price of employing 
Poisson brackets with the volume functional. 

Let us see this in more detail^^: 

We will consider the simplest case: The isotropic, fiat FRW cosmologies for which the line element 
reads 

ds^ = -dt^ + a{tfdx^ (2.16) 

where a is the scale factor. Here we take t, a to have dimensions of length. One fixes both the spatial 
diffeomorphism constraint and the Gauss constraint by setting 

4:=g5^/3, El.= p5'; (2.17) 

where p, q are spatial constants. It follows that a = \/\p\ and that {p, q} = 8tc(3G are canonically 
conjugate. It is easy to see that with ()2.17|) both the Gauss and spatial diffeomorphism constraint 
vanish identically as they should and that the gravitational contribution to the Hamilonian constraint 
becomes C = —aq'^^/\p\ with a positive constant a. For a (homogenous) scalar field we find 

^scalar ^ ^ /H V(0) (2.18) 

As one can show (see the end of section 0} the classical singularity is encoded in the kinetic term 
only. 

The holonomy and flux become respectively 

A{e) = cos(gre) + r.niMreq), Ef{S) = [ [ ds^Hp (2.19) 

Js 

where e"(l) — e"(0) =: Tgn^, (^i)^ = 1- By varying e and taking traces we may generate the 
functions := exp(irg) where r can be any real number. The reduced holonomy flux algebra 
therefore is generated by the and p. 

The reduced Hilbert space is now obtained by asking that the form an orthonormal basis, i.e. 
< Tr,Tr' >= 6r,r' SO that they play the analogon of the SNF. Here the continuous graph label and 
the discrete spin label of the SNF have joined into one continuous label r. The operators T^, p act by 
multiplication and derivation respectively, explicitly T^/ = T^+r' and p := ih{p,Tr} = —ripPTr 
where = SnGh is the Planck length squared. 

^^Originally one tried to define LQG as a symmetric sector of distributions on a dense sunspace of the LQG Hilbert 
space together with the dual action of the full LQG operators. However, this dual action does not preserve the 
symmetric sector. See also the end of this section. We adopt here the correct and modern point of view spelled out 
in 137]. 

""^^What follows is a drastically simplified version of |37j . 
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Finally, in order to quantize objects like the inverse scale factor 1/a = we notice the 

classical identity 

1 T,JT_,„ VM} - T,,JT,„, /H} (2.20) 



\p\ iSnpGro 



which holds for every ro 7^ 0. Notice that ^/\p\ oc V^^^ where V is the volume^^ so that ()2.20j) is 
indeed analogous to ()2.14|) . Notice that the spectrum of the volume operator is essentially given 
by |rp/^ with Tj, as normalizable eigenf unctions (it takes a continuous range but consists only of 
eigenvalues, there is no continuous part). Hence the only zero volume state is Tq = 1. It follows 
immediately that if the inverse scale factor can be densly defined on the Tg then at zero volume it is 
bounded. 

Expression ()2.20p is readily quantized by substituting Poisson brackets by commutators divided 
by ih. This gives 



J. ^ Tr, ^/\p\T-ro - T_ro \/\p\Tro ^ ^/\r + Tq] -y\r-ro\ ^ 



which is bounded (in fact vanishes) at r = 0. It is even bounded for all r and coincides with 
the operator for r ^ tq. The parameter tq is a quantization ambiguity, called a factor 

ordering ambiguity which is also present in LQG |SH| although it is there labelled by a discrete 
rather than continuous parameter. Notice that by the constraint C = Cgeo + C scalar = we obtain 
oc p^/(2a^) + Vo? which means that q cx at a — >• (see the discussion at the end of section Hj) 
so that Cgeo 

oc a classically. However, quantum mechanically Cg^o is bounded at a — > because q 
is quantized by the bounded expression (T^q — T^ro)/ (2«ro). 

So far we have only discussed the isotropic model. In the non - isotropic but diagonal model [201 
similar remarks apply because the three scale factors mutually commute. Hence the inverse scale 
factors can be treated separately just as in the isotropic case. 

Let us summarize: 

LQC is the usual cosmological minisuperspace phase space quantized by LQG methods. It borrows 
in an essential way the Hilbert space |HI and the key technique [T2j. LQC therefore seems to to 
confirm the aspects [3112] for LQG. 

However, despite these similarities there are important differences: 
To begin with, LQC does not have a graph label. That is of course expected because in cosmology 
all points are equivalent due to homogeneity, hence there is only one point and thus no room for 
a graph. Consequently, there is no such thing as the valence of a vertex. These facts lead to an 
enormous simplification of the volume spectrum which is only modestly degenerate (multiplicity two 
for all r 7^ 0, multiplicity one for r = 0) while in LQG the volume spectrum is infinitely degenerate, 
in particular for zero eigenvalues (just apply a diffeomorphism to a given volume eigenstate which 
changes the graph). Moreover, the volume spectrum in LQG has discrete range while in LQC it has 
continuous range. Next, in LQG the Hamiltonian constraint modifies the graph of a spin network 
state while in LQC that is not possible since there is no graph. That leads to an enormous simpli- 
fication concerning the computation of the kernel of the Hamiltonian constraint which boils down 
to the solution of a recursion relation within a finite dimensional parameter space while in LQG the 
corresponding parameter space is infinite dimensional. Finally, while it is possible to interpret 

^•^For fc = this volume is actually infinite. We mean here the volume divided by the comoving volume. 
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the LQC states as distributional LQG states [HH] (simply multiply an ordinary LQG state by a 5 dis- 
tribution supported on homogeneous connections) this distributional space of states is not preserved 
by the (dual) action of the operators of LQG, hence LQC is not an invariant (distributional) sector 
of LQG. It really is a model whose connection to LQG involves the classical truncation, thus inho- 
mogeneous fluctuations are switched off by hand rather than being quantum mechanically suppressed. 

We conclude: By definition of a toy model with a finite dimensional phase space, LQC cannot 
model the QFT aspects of LQG. LQC models very successfully two aspects of LQG, namely those of 
[HlEl- However, it is far from granted that those predictions of LQC which are sensitive to the QFT 
aspects will hold in LQG as well. As we just saw, the spectrum of the inverse scale factor is sensitive 
to those aspects since the volume operator is. It is precisely one of the purposes of the present paper 
to investigate whether these differences lead to a qualitatively different picture in LQG when applied 
to cosmology. 

3 Unboundedness of the Analogon of the Inverse Scale Fac- 
tor in Full LQG 

We can explain the reason for the unboundedness of the operator concerning the inverse scale factor 
in very intuitive terms without going through the (admittedly lengthy) analysis. We will compare 
this to the isotropic LQC model, similar arguments apply to the the diagonal models: 
The spectrum of the inverse scale factor can, roughly speaking, be obtained by taking the difference 
between neighbouring eigenvalues of the volume operator (discrete derivative). This is the case 
because, according to ()2.14|) the inverse scale factor can be obtained from the Ko - Dreibein as 
e = {A, V} where V is the volume and A is the connection. This is the key technique of jT3] as 
mentioned above. Upon quantization the connection gets replaced by a holonomy, that is, e oc 
h[h~^,V] = V — hVh~^ where /i is a the holonomy of A along an edge. The holonomy changes 
the spin labels of a spin network state by one unit, hence we get a discrete derivative in the spin 
parameter space. Now the spectrum of the volume operator depends parametrically, among other 
things, precisely on the spins of the edges of a spin network attached to a vertex and additional 
intertwiner quantum numbers and this is the discrete derivative of the volume spectrum that we 
mentioned. 

If we view the spectrum of the volume operator as a function of these parameters then in isotropic 
LQC this parameter space is one dimensional, it is a function j ^-^ (in LQC the parameter 

j := r of the previous section takes continuous and not discrete values but that does not change the 
qualitative picture). In LQG on the other hand it is arbitrarily high dimensional, the dimensionality 
being governed by the number n of vertices and the valences of the individual vertices, it is a 
function of the form ••, jiT"*) ^ S"=i ^viJi'K jj^'*) where the n, iV„ can be arbitrarily large. In 
both cases the functions V are nowhere singular except when at least one of the j tends to infinity. 

The point is now that in isotropic LQC the zero volume condition V{j) = has a finite number 
of solutions j = jo < oo (in fact one). The spectrum of the inverse scale factor at zero volume is 
then, roughly speaking, given by V{jo + 6j) — V{jo) (using suitable discrete units for the parameters 
j, in LQC 6j = tq) which is finite since jo < oo whence we get boundedness at zero volume. In 
LQG on the other hand the situation is more interesting: The zero volume condition now gives an 
infinite number (parameter space of codimension one) of solutions. Pictorially speaking, the spectrum 
of the volume operator plotted against the parameter space can be envisioned as a complicated 
landscape of mountains which is disrupted by many valleys of zero volume which intersect each other 
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in a complicated way. The boundaries of those valleys have arbitrarily steep inclination and the 
mountains are arbitrarily high. This is because the solutions of the zero volume condition are now of 
the form = ■■iTmI) ^"^^ •^2'''*? --^jNl ~^ 00 in many ways (along many valleys). 

It is clear that now for each v V;((jf'' = jo\j2\ --JnI)) - ^viiJi^ = jo,j2\ ••, j£) + Sj^"'') with 
Sjl^^ = 6ki, I E {I, .., Ny} can become arbitrarily large. This is what we are going to prove explicitly 
in our companion paper [28j based on the recent simplification of the analytical expression for the 
matrix elements of the volume operator 

One may object that the diverging configurations {ji'\j2'\ jj^'') for each individual v do not 
look isotropic in general and that one should require isotropy as in LQC which perhaps should enforce 
j^^ = .. = which could lead to a bound since the parameter space is now one dimensional^^. 
However, here we now invoke a second key difference between LQC and LQG: LQC is a quantum 
mechanical toy model of LQG. It can be thought of as LQG where all the vertices of a spin network 
state are to be thought of as equivalent due to homogeneity. In other words, LQC is like LQG except 
that we only consider spin network states with a single vertex (and a few edges, usually no more 
than three). LQG on the other hand is a field theory and spin network states can have an arbitrary 
number of vertices. Consider now a zero volume spin network state Tx^^ where A is the eigenvalue 
of the inverse scale factor and f is a vertex at which the graph 7^ with a single non co - planar^^ 
vertex v of the spin network state is located. Let 5* be a finite^^ set of n vertices obtained from a 
random sprinkling of points into the spatial manifold PU] and consider the state ^p := ^vesT\^v We 
may arrange that the 7„ are mutually disjoint and that their orientation is random as well. This 
state is homogeneous and isotropic on large scales but not on microscopic ones (the violation scale 
of homogeneity and isotropy depends on the number of points in 5* and decreases with n). The state 
is normalized and still an eigenvalue A eigenstate of the inverse scale factor. This proves that we 
may construct kinematical states with large scale homogeneity and isotropy in LQG with arbitrarily 
large norm of the inverse scale f actor 

Let us now go into more details: 

In order to match the LQC calculations as closely as possible, what we really have studied is the 
contribution of a scalar (or infiaton) field to the Hamiltonian constraint (smeared energy density) 
operator whose classical expression is 



2CscaiariN) = d'xN [^=== + ^dct (g) (g'^V.a^.fc + 2V(0))] (3.1) 

Here is the scalar field, vr its conjugate momentum, g^fe is the pull back to the spatial manifold a 
of the four metric g^u on the four manifold M = R x a and is a test function (see also section |2)). 
In an exactly homogeneous situation, say the FRW metric, the derivative term drops out and the 
potential term becomes proportional to while the kinetic term becomes proportional to where 



^"^For the very special case of matrix elements we are going to use a similar expression was also derived in |36| using 
a different method. 

^^Indeed, if we choose ji = j2 = = j in formula (|3.5|l then the result is bounded in j. 

^^The volume operator vanishes on (non) gauge invariant vertices which are not at least (tri) four - valent and at 
which not at least three tangents of the adjacent edges are linearly independent. 

^^This will be sufficient if the spatial manifold is compact as ioi k — 1 FRW models or if we are interested only a 
compact, observationally accessible (inside the Hubble radius) region. In the fc = 0, ~1 case we can use in the infinite 
tensor product |1T] . 

-•^^This is a purely kinematical result as in LQC. In fact it must be because the inverse scale factor itself is not a 
Dirac observable, see next section. 
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a is the scale factor of the FRW metric. The classical initial singularity limt^o'^('^) = is therefore 
prominent in the kinetic term. Of course in full LQG the kinetic term must be and can be considered 
as well, see [12], but it will be sufficient to consider the kinetic term to demonstrate our point. 
Consider therefore the kinetic term 

'^C^ZariN) = I d'xN (3.2) 

Its quantization is derived in [T3] on the kinematical Hilbert space Ttkin = 'Hgeo ® 'Hscaiar- An 
orthonormal basis in this Hilbert space jl2] consists of the charge - spin network states T^^ f®Tv{'y).fi- 
Here 7 is a graph with vertex set ^(7), j denotes a labelling of each edge e of 7 by half integral 
spin quantum numbers je, and /, n respectively denote a labelling of each vertex f of 7 with a 
gauge invariant intertwiner 1^ and integer respectively. The action of the corresponding operator, 
densely defined in the finite linear span of this basis, can now schematically be written as 

2^.(iV)7;j/,.^= ^(^)[^^M^.J/] ® [^(^)'^n7),n] (3.3) 

^ey(7) vclet^gj 

The explicit meaning and action of the genuine operators (not operator valued distributions) can be 
found in 01] and is also reviewed in [2H] for the convenience of the reader. For the purpose of 
this paper it is enough to report the result of our calculation. 

The analogon of the LQC inverse scale factor (cubed) in LQG is of course the non - negative 

operator , ^ iv). To show that it is unbounded from above on zero volume eigenstates it is 

sufficient to construct one single system of states whose members have arbitrarily large norm. The 
simplest zero volume eigenstates are spin network states 7 f which have at most three valent 
vertices. They are very simple because the intertwiners J„ are then uniquely determined. We are 
thus interested in the numbers 

T.:^\?=<T,.lA^^{v)fT^3> (3-4) 



where w is a three valent vertex of 7 whose adjacent, non co - planar edges carry the spin quantum 
numbers ji < J2 < js subject to the constraint that J3 G {ji + ^25^1 + 32 — 1, ••5 lii — i2\\- To make 
things definite, consider a graph 7 with three edges and two three valent, non co - planar vertices, 
see figure 01 




Figure 1: Simplest 3 va- 

lent gauge invariant graph 7: 

1/(7) = {vi,V2], E{-i) = {61,62,63} 



15 



Remarkably the number (j3.4p can be computed analytically. This is possible due to a recent 
improvement in the technology for the computation of the matrix elements of the volume operator 
for full LQG j^. The analytical formula for ()3.4|) is given by the surprisingly compact expression 



C2 



X 



X 



32 



9(l + 2ii)'(l + 2j2)'(l + 2j3)' 



-3 (2 (-l)^^'^ + {-if'-^yA^ {A, -A, + Asf 
-3 (2 (-1)^^'^ + (-l)^^'^) ' A, {-A, + A, + A, 
-3 (1 + 2 i-lf^'^^^'^y As (A + A,- Asf 



2(_i)2i2 + (_l) 
-Ai + A2 + As) {Ai -A2 + A-i) (ii + A2- A3) ] 



2i3 



(3.5) 



where C is a positive numerical constant depending on the regularization of the operator 



and the volume operator. Furthermore (ji, j2, js) are the spins of the 3 edges (61,62,63) outgoing at 
the vertex v and Ak := JkUk + !)• Moreover 

1 

Via = [(-ji + j2 + h + 1) (ji - 32 + js) (ji + J2 - js) (ji + 32 + h + l)Y 

VlB = [(-jl + J2 + Js) (jl - 32 + J3 + 1) (jl + 32 - h + 1) (jl + 32 + h + 2) 
V2A = [(-jl + 32 + J3) (jl - j2 + 33 + 1) (jl + j2 - js) (jl + j2 + js + 1) 
V2B = [(-jl + j2 + js) (jl - j2 + js) (jl + j2 - js + 1) (jl + j2 + js + 2) 
VsA = [(-jl + j2 + js) (jl - 32 + js) (jl + j2 - js + 1) (jl + j2 + js + 1) 
VsB = [Hi + 32 + 33 + 1) (ji - 32 + 33 + 1) (jl + 32 - Js) (jl + 32 + 33 + 2)J ' 
We now display a few diverging series of configurations (ji,j2,js) where we always plot 

^, neglecting the pref actor 
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"Oscillating" 



Jl = J2 



11 
2 



If we set jl = j2 = f- where e N, we get: 



Vdet((?) 



128 V2 (-1 + (-1)^^) js' (-3 (2 + (-1)-'^) + (1 + (-1)'^'^) is) 



9 (I+J3) (l + 2j3)' 




Asymtotically this increases as 

II '^~l v)T V - 128.6.V2.j3^ 



OC 



9(l+j3){l+2j3)^ 

''oc 35.9 -ia^ 



Figure 2: Plot for ji — j2 — ^ where 
e N with 1 < is < 40. The graph 
oscillates between (if js even) and an 
increasing value (if js odd) 



(3.6) 



Jl - J2 - Y + 2 



If we set jl = j2 = Y + I where ja e N, we get: 



X 128V2 (l + (-ir) ,73'(l+i3)' (-3 (-7 + 3(-l)^^ + (-lf^^) + (-l + (-lf^^) is) 

x(- (2 J (is^ (l+i3))^) + ((l+i3)' (3 + 2j3))') 

(j3(l+j3)')' -(j3 (3 + 5j3 + 2j3'))' 




Asymtotically this increases as 



v)T^j\\ (X 5.9-10-^ia^ 



Figure 3: Plot for ji = j2 = where 
ja G N with 1 < ja < 40. The graph 
oscillates between (if ja odd) and an in- 
creasing value (if js even) 



(3.7) 
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Increasing 



ji = I J 2 = ja + I 



If we set ji = |, j2 = is + ^ where js G N, we get: 



Vdct(g) 



i\ 2 



3 (I+J3) (1 + 2J3)' 

x4i3 (-9 + 21j3 + 14j3') (-(35O3 (2 + j3))^) + (-3 + 4^3 + 4^3') ^) 
X (-2 (i3 (2 + ^3))^ + \/2 35 (-3 + 4j3 + 4j3')^)' (O3 (3 + 2j3))^ - (-6 + 9^3 + 6^3')*) 



IN 2 



, ^ (v) T 




Asymtotically this increases as 



•^3 II 1 f T -I 



oc 1.06- 10"^ jg^ 



Figure 4: Plot for ji = |, j2 = js+l where 
js G N with 1 < js < 40. The graph first 
decreases for 1 < js < 3 and is for js = 3 
. It increases for js > 3 



Identical 

If we set ji = j2 = js and more general ji, j2, js G N (all spins integer numbers) in ()3.5|) then 
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General Configurations 

Finally we plot the numbers ()3.4p as ji,j2 vary over their allowed values at given js for various values 
of ja. The landscape of low valleys together with their arbitrarily high and steep boundaries is nicely 
illustrated: 

Taking the general result ()3.5|) (without the pref actors (£p)^|Z|t) we use the quantity 



Q 



[l I ^dct(g) ^'"'^ ^ijllO'i'ia, is)] |ji -j2\ < h < ji +32 and ji +j2+h is integer 

1 1 -7^77=^ (^) ^7jl I (ii ' ^2, J3) ^ 



y/det{q) 



else 



where the 30 is arbitrarily added in order to shift the plot upwards. We make a three dimensional 
plot (in the range ^ < ji < jmax, | < J 2 < jmax for each fixed value 5 < J3 < ^ : 

I turns out that the non-zero configurations are grouped symetrically along lines parallel to the 

ji = 32 axis. The reason for this is of course the integer requirement ji + j2 + is = integer. Therefore 
we will get contributions on the ji = j2-axis only if is integer. Because (13. 5p is symmetric with 
respect to the interchange of ji ^ j2 we may restrict ourselves to the range ji > j2 ■ 

Additionally we contribute for each of those 3D-plot a 2 dimensional plot along the lines 

1 1 

j2 = ji-l mit < / < mm[j3,jrnax - 2^ ' 2--^^- ^""""^ ^ 

The restriction for the parameter / is a result from the requirements |ji — j2\ < Js < Ji + J2 from 
which for ji > j2 we may remove the modulus. In order to give a better impression we have joined 
only non- vanishing values of Q along the lines described above^^. 



^^For integer J3 also every second configuration on the lines gives a (the identically case). Hence, if we would 
join all the points, the plots would oscillate between and the now plotted curves and it would be hardly possible to 
see anything useful from them. 
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5 10 15 20 25 



Figure 9: Plot for js — 6 Figure 10: Plot for — 6 




Figure 11: Plot for = ^ Figure 12: Plot for js = y, j2 = ji - I. 

The different curves are (bottom to top) : I — 

13 13 
2' 2' ■ ■ ■ ' 2 




Figure 13: Plot for = 7 Figure 14: Plot for ja = 7 
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Figure 15: Plot for J3 = ^ 



Q 




Jl 



Figure 16: Plot for = ja = ji - I- 
The different curves are (bottom to top): I = 

13 15 
2' 2' • • • ' 2 



Due to the large kernel of the volume operator even for four - and higher valent vertices jHSj a 
qualitatively similar behaviour will prevail in those cases as well. The unboundedness of the "inverse 
scale factor" operator in full LQG, even on homogeneous and isotropic states is herewith established: 
Let 7i, ..,7„ be a random distribution of mutually disjoint but diffeomorphic graphs, say of the type 
of figure 01 which is homogeneous and isotropic with respect to the comoving spatial metric of a 
FRW spacetime on a scale large compared to the Planck length within a compact region (the smaller 
that scale the larger n). The combined graph 7 := U^^^^ jk looks then homogeneous and isotropic. 
Consider the state of unit norm 

T,,A := ®Li T^,J (3.8) 

where we choose the same triple of spins for every graph and such that the numbers (j3.4j) are equal 
to the same value for every vertex of the 7^. Then obviously all the norms (|3.4p with T^j chosen 
as the homogeneous state ()3.8p equal A^ for every of the 2n vertices of 7. But A can be arbitrarily large. 



Let us compare this with the situation for the hydrogen atom: In both LQG and LQC the spectrum 
of the inverse scale factor is purely discrete, so we should compare with the bound states of the hy- 
drogen atom. Restricted to the subspace of bound states it is indeed the case that the Hamiltonian 
of the hydrogen atom is bounded from below and above and this is not changed when we switch from 
the finite number of degrees of freedom of QM to the infinite number of degrees of freedom of QED. 
However, the interplay between LQC and LQG is different: In LQC we have boundedness, but after 
switching on LQG we get unboundedness. 

4 Boundedness of the Inverse Scale Factor Expectation Val- 
ues at the Initial Singularity 

The result of the previous section sounds rather negative at first. However, it is not as we will 
now explain. It just shows that the boundedness of the inverse scale factor cannot be a 
mechanism for avoidance of the local singularity in LQG. 

Namely, what we have shown is unboundedness on certain homogeneous and isotropic spin net- 
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work states. However, that does not mean that the inverse scale factor is unbounded on all homoge- 
neous and isotropic states. In fact, we should worry only about boundedness on those states which 
describe a collapsing universe which is homogeneous and isotropic on large scales. Indeed, in this 
section we will display the result of a calculation which indicates that the expectation value of the 
inverse scale factor with respect to (kinematical) coherent states flB^ peaked on any homogeneous 
and isotropic point in phase space (e.g. a FRW universe along the classical singular trajectory) is 
bounded from above in LQG even at the Big Bang. 

The details of our calculation can be read in the companion paper j2H]- However, again we can 
explain the basic mechanism for this boundedness result in very intuitive terms: 
Spin network states are, roughly speaking, eigenstates for the kinematical geometric operators such 
as the triad operator out of which the inverse scale factor is assembled. They are thus maximally 
sharp for half of the degrees of freedom (the electric fields) but maximally unsharp for the other 
half (the conections) which are encoded in the holonomy operators. Coherent states intermediate 
in between those two uncertainties and in fact minimize the Heisenberg uncertainty bound as far 
as conjugate canonical pairs are concerned. This is achieved by superimposing basis states such 
as spin network states with carefully chosen coefficients. These coefficients depend on a point in 
the classical phase space, in our case a field configuration {A, E) of connections and electric fields 
respectively. The coherent states are thus labelled by these points and they are peaked in phase 
space in the sense that the overlap function | < iIj{a,e)i'^{A' ,e') > P vanishes exponentially fast unless 
||y4 — A'lp, — < s where ||.|| is a suitable norm on the classical field space and s is a 
parameter which encodes the Gaussian width of the wave packet or the volume of a cell in phase 
space (it is proportional to h). The Gaussian decay of the overlap function implies that in terms of 
spin network states the state is of the form ip{A,E) = J2j (^5 ^) where we have not displayed 
the graph label for simplicity and where the Cj decay (symbolically) as exp(— sj^). The label j is a 
compound label for all participating spins and intertwiners. Now while the norm of the inverse scale 
factor in the states Tj is unbounded from above, that bound only diverges with a finite power as 
j — >■ 00. This divergence therefore is completely swamped by the Gaussian decay of the exponential. 
This holds for rather generic values oi A,E and in particular when we choose a homogeneous and 
isotropic point in the phase space. Hence the mechanism is similar as in statistical physics: "Energy" 
(the argument of the Gaussian factor) wins over "entropy" (the number of configurations leading to 
divergent eigenvalues). This a rather different mechanism which does not require boundedness of the 
inverse scale factor. 

Let us go into more detail: 

Given a point (Aq, Eq) in the classical phase space, our (not normalized) coherent states are of the 
general form 

^(AcBo) = X] ^7,(Ao,i?o) := [e~*^^^ ^A']a'^a«(Aq,Eo) (4-1) 

7 

where the uncountable sum extends over all finite graphs and ip'y^{Ao,Eo) is a countably infinite coher- 
ent superposition of spin network states over 7. The coefficients in that superposition are uniqely 
determined by the generator C of the coherent states, also called the complexifier |17j, which is 
supposed to have the dimension of an action. In fact, as we can see from ()4.H) . the coherent states 
are nothing else than the heat kernel (with respect to the operator C) evolution of the 5 distribution 
with support at A which is then analytically extended to A^{Aq,Eq). The complexified connection 
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also is generated by C via 

^= E ^ {C, = A + z {C, A} - - {C, {C, A}} - '-{C, {C, {C, A}}} + ... (4.2) 

n=0 

This may look unfamiliar, however, this is precisely the way that the unnormalized coherent states 
for the harmonic oscillator are constructed for which the complexifier reads C = p'^/{2muj) and 
the classicality parameter is s = h/ {mui) (the role of the connection is of course played by the 
configuration coordinate q). 

Ideally, according to the subsequent section, C should be adapted to a given Hamiltonian, but 
as we said, here we will carry out a completely kinematical calculation to get a first insight. It is 
therefore appropriate to choose C in such a way as to obtain the simplest possible states. One such 
choice involves the square of the area operator |37j which in turn is proportional to E'^ and thus 
will give rise to states very similar to the ones for the harmonic oscillator because in both cases the 
complexifiers are of the momentum squared type giving rise to Gaussian wave packets. Rather than 
defining C we directly display the corresponding states which in that case are, roughly, of the direct 
product form 

oo 

^7,(Ao,i^o) = n ^-'(^.^) ^^^"^ ^e,(Ao,Eo)(^) = + 1) e-^^'(^'+i) Tr(7r,(^/e(A, E) A(e)-i)) 

eeE(7) 2j=l 

(4.3) 

Here vTj is the spin j irreducible representation of SU{2). Notice that that the wave function 
'4'e,{Ao,Eo){A) is a function of A (or of the SU{2) valued holonomies A{e)) which depends para- 
metrically on {Ao,Eo). The dimensionless parameter t{e) depends on the details of C and satisfies 
t(e) > 0, t(eoe') = t(e) +t(e'), t(e"^) = t{e). Finally, ^e(^o,^o) ~ exp(zrjE^(5e)/2) Ao(e) where Se 
is, roughly speaking, the face of a cell complex dual to 7 (i.e. each intersects 7 only in an interior 
point of e) which also depends on the details of C). Experts will see the representation theory of 
SU{2) and the Peter&Weyl theorem at work in (j4.3|) . 

As explicitly shown in expectation value computations with coherent states are qualitatively 
unaffected when we replace the group SU{2) by f/(l)^ for which however the analysis simplifies 
tremendously. The corresponding coherent states are then given by 

00 3 
^e,iAoMA)= E e-'^^Ur^ ll(9iiAE)A^{e)-'r (4.4) 

|niNri2|,|ri3|=l j=l 

where A^{e) = exp(i J^A^) and where 

giiAo, Eo) = exp( / dS^E^^ + i ! A^^) (4.5) 

To carry out concrete calculations we just have to insert concrete values for Aq, £"0 into those formulae. 
In cosmological applications we choose, as in LQC (see section |21), e.g. Aq^ = q5{ and E^^ = p5{ 
where a = \/\p\ is the scale factor. Thus in LQG, rather than performing a cancellation of the 
inhomogeneous degrees of freedom as in LQC, we perform a suppression of their fluctuations by 
using appropriate states. 

What one actually does in the computations is not to use the states ip(Ao,Eo) but rather the "cut 
- off" states ip'y^(Ao,Eo) foi' finite but sufficiently fine 7. Roughly speaking, the analysis of ^\ ITT] 
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shows that the optimal graph, for given resolution scale L of a measurement^^ should have an 
average edge length scale (with respect to the metric to be approximated) of the order of a geometric 
mean e ~ L^tp^ where < A; < 1 is a parameter that depends on the (partial) observables to 
be approximated. The reason for working with the cut - off states is that the iP{Aq,Eo) are not 
normalizable on the kinematical Hilbert space due to the sum over the uncountable number of 
graphs. They merely serve as a generator for the cut - off states which are normalizable. 

What we have done then is to estimate from above the expectation value not only of the inverse 
scale factor in LQG for a general valence M of the vertex for arbitrary (Aq, Eq) but also the most 
general polynomials built out of generalized co - triad operators coming from classical expressions of 
the form {{A, V^})^ which become upon quantization operators of the form [-^j Y\j=i Q.iji'^-:^) = 
[^]^ ri/Li ^e_r [(^e/)~^ ^""l where r, N are fixed, positive, rational and natural numbers respectively 
which are determined by the concrete coupling between matter and gravity Here h{ are the 

holonomies of the connection A along edges of the graph, the classical volume expression V turns 
into the volume operator V and Poisson brackets into commutators. A bound is given by 
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where M denotes the valence of the vertex, the degree of the polynomial, Z is a constant 
numerical prefactor, depending on the Immirzi parameter and the regularization of the volume op- 
erator. 

Moreover 

min{T/} A 



T :-- 
p : = 



1 + f 



max{|p^(m)|} 



with Tj = ^Jt{el), t{ei) is the classicahty parameter mentioned above andpj(m) ^ -p- dSaEj. 
Here m = {Aq, Eq) = {A, E) is the point in classical phase space the coherent state is peaked on and 
L, e respectively are the resolution scale and graph scale mentioned above. 

Specializing to homogeneous data m = {A^^, E^-) = {qSi,p6j) gives p = p ■ maxCj'^""**'^"'^*'*, with 

numerical constants (^j™"''*''"^*) carrying information about the details of the coherent state. 

Taking the limit a — > which causes p = 0, A = 1 demonstrates that nothing dramatic happens 
and we get the huge but finite bound 
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(4.7) 



For the case of the inverse scale factor we have to set (r = 1/2, = 6) and to multiply by 

1 _ 1 



{£p)12- 



^"^Given a metric go to be approximated, we want fluctuations in the length to be much smaller than L where L is 
to be measured with respect to go. On the other hand, the Planck scale discreteness of LQG forces e ^ L. 
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Neglecting all terms of order t° 
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These results prove our claim: The expectation value (and thus also the norm, just double N) 
of the inverse scale factor is bounded with respect to kinematical semiclassical states even when the 
classical metric that labels the coherent state becomes maximally degenerate. As in LQC this is a 
strong indication for a drastic modification of the effective Friedmann equations near a = \2'2\ . 
This should be worked out in more detail in order to decide, along the scheme given in section El 
whether the initial local singularity is avoided or not. 

Remarks: 
1. 

We have computed the expectation value of the analogon of the inverse scale factor with respect to 
kinematical coherent states whose peak follows the classcial (singular) trajectory. One should do the 
same with respect to the matter degrees of freedom. However, it is clear from the outset that this 
will not change our result: The classical Friedmann equations in the gauge outlined above are 

- '^^''^ + ^) + A - |. * ^ -a^V (4.8) 

where if = - is the Hubble function and vr = a^<p the conjugate momentum. Consider the case of 
a potential which depends polynomially on but not exponentially. We want to solve ()4.8|) close 
to a = 0. We now make the self-consistent assumption (to be checked) that a^V (and thus a^V) 
vanishes in the limit a —* 0. Hence vr becomes constant close to a = according to the second 
equation in (|4.8p . Then we may neglect the potential {V), cosmological (A) and curvature (k) term 
in the first equation as a — and thus we may solve a{t) oc t^^^ and 0(t) oc ln(t). Hence indeed 
a^(f)"- oc t[ln(t)]" vanishes at t = and our assumption was correct. It follows that the expectation 
value with respect to coherent states, taking into account both matter and geometry, of the scalar 
contribution to the Hamiltonian constraint ()2.7|) . whose homogeneous reduction is 7r^/(2a^) +a^V, is 
indeed non-singular even when the coherent states are evaluated along the exact classical trajectory, 
if and only if the expectation value of the analogon in ()3.3j) of is bounded from above. This is 
precisly what we verified. 
2. 

We have shown unboundedness of the inverse scale factor in LQG (with gauge group SU{2)) and 
did expectation value calculations therefof for an approximation to LQG using U{1)^ instead of 
SU{2). It is maybe worthwhile to mention that the inverse scale factor for that approximation is 
also unbounded on zero volume eigenstates although one has to use graphs of valence at least five 
here which is due to the detailed difference between the non - Abelean group SU (2) and the Abelean 
group f/(l)^. 
3. 

Notice that the bound in formula (j4.6|) is completely general. It does not depend on the value of 
the classical connection A that labels the coherent state, it only depends on the electric field E. 
Setting = as appropriate for the isotropic singularity gives the upper bound ()4.7p . However, 
we may consider any other homogeneous model, for instance inhomogeneous ones such as Bianchi I 
(Kasner), Bianchi IX or Mixmaster Models (BKL scenario). In these cases the different components 
of the three metric vanish with different rates as we approch the singularity (or do not vanish at 
all or might even diverge, at least in the absence of matter). As we can see from the above general 
estimate, the important question is whether the quantity p defined above, which is proportional 
to = 7det(g)e», remains finite in all those cases. For instance, in the Kasner models with line 
element ds'^ = —dt"^ + Yyl=i t'^^°{dx"'y where the real parameters Pa are subject to ^^Pa = ^aPa — ^ 
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we find in the diagonal gauge Ej = 6jt^~P"- which stays bounded for t ^ because Pa < 1- In 
particular, contracting directions (with increasing t) are allowed. This covers all Kasner models and 
shows that the expectation value calculation performed here also hints at local singularity avoidance 
for at least some class of anisotropic models. Whether this holds for all homogeneous Bianchi models 
and also after matter is coupled is a task left to the reader. 

5 Scheme for Determining the Presence or Absence of the 
Initial Singularity 

The computation of the previous section is very encouraging. Unfortunately, it does not yet prove that 
in LQG the Big Bang singularity or any other classical spacetime curvature singularity is quantum 
mechanically absent. As we have explained at length in the introduction, fundamentally, singularity 
avoidance must be formulated in terms of physical states^^ and physical operators which has not 
been done so far even in LQC^^ because the volume, the inverse scale factor and the value of the 
infiaton field etc. are not gauge invariant quantities (they do not commute with e.g. the Hamiltonian 
constraint). 

Fundamentally, what one would need to do in order to establish a singularity avoidance result is 
something along the following lines^^: 

1. Ideally one should have sufficiently explicit knowledge about the physical Hilbert space Tiphys 
and physical states. An avenue towards this goal is the Master Constraint Programme |17j . 
Notice that physical states automatically solve the Quantum Einstein Equations. That means 
that the quantum gauge evolution ^Hl, defined similarly as in as the coefficients of the 
physical states when expanded into (spatially diffeomorphism invariant) spin network functions, 
is non - singular by definition thanks to the fact that operators like the inverse scale factor 
are densely defined in LQG. Moreover, in LQG there are physical states which are composed 
of spin network states that describe a sign fiip with respect to the expectation value of the 
triad orientation sign operator 03] and thus may describe quantum geometries which include a 
Pre Big Bang in the sense of ^H]- As we have already outlined in the introduction, the global 
singularity is avoided if and only if there are sufficiently many semiclassical of these physical 
states. 

2. (In)famously, canonical quantum gravity, especially in the cosmological context, is a theory 
without a canonical Hamiltonian. However, one can construct physical Hamiltonians using 
the relational point of view (23 EO] . These are Dirac observables by construction. They are 
obtained by i) using non - gauge invariant quantities P, T where P is called a partial observable 
and T a collection of^^ clock variables, ii) computing the gauge fiows (Hamiltonian evolution 
with respect to the Hamiltonian constraints^^) t ^ at{P), at{T), iii) choosing parameters 

^^The singularity effects A and B for the global singularity have indeed been discussed for solutions of the Wheeler 
DeWitt equation in LQC. However, that does not grant that a sufficient number of them are normalizable or have non 
vanishing norm in the physical inner product. Once again, formal solutions are not automatically physical states. 

^•^In terms of the the partial and complete observables discussed below, there is in general no relation between 
the spectra of the complete and partial observables Hence kinematical boundedness does not imply physical 

boundedness as far as the local singularity is concerned. 

■^''This is just a sketch of ideas, not a "paradigm" . 

^^In general one needs as many clock variables as there are constraints. 

^^There are also as many unphysical time parameters t as there are constraints. In GR they become functions and 
are called lapse and shift. 
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r and solving the equations at{T) = r for t and iv) inserting those values of t into at(P). 
The resulting object Pt{t) has the physical interpretation of giving the value of P at the 
moment (gauge parameter t) when T has attained the value r. The function Pt{'t) is called the 
corresponding complete or Dirac observable. The physical Hamiltonians Ht{t) are in general 
explicitly dependent on the physical time r and are defined by {Ht{t) , Prij)} := d/drPTij) 
for all P. Of course there are many possible choices for Ht depending on the choice of the 
"clocks" T (in GR we need an infinite number of those). In cosmology it would be nice to 
use the scale factor (total volume of space) as a clock. If a corresponding positive (or at least 
semibounded) and physical time independent Hamiltonian exists and can be defined as a self 
- adjoint operator on the physical Hilbert space then we may use it to define evolution on the 
physical Hilbert space. 

3. If a physical Hamitonian is given, then we will need coherent states in Tiphys which are preserved 
by Ht for sufficiently long physical time. A general idea for how to do that has been given by 
Klauder at least for systems with a finite number of degrees of freedom. 

4. We are now in the position to address the issue about the local singularity: Given a physical 
coherent state ip peaked, on scales large compared to the Planck scale, corresponding to r = 
tq ^ Tp, on classical singular initial data, consider the map r i-^< ip,P'r{T — tq)iP >phys=< 

V^(t— To), Pt(to)^/'(t— To) >phys with iP{t) := e^'^^^ip. Let P be any classically singular quantity, 
say the matter energy density or a curvature scalar density (suitably smeared). Then the local 

singularity, with respect to P, is absent provided that < ip , Pt{—tq)%Ij >phys is finite (assuming 
that r = corresponds to vanishing classical volume which will be the case if we choose V as 
one of the clock variables). 

Notice that this ambitious scheme has not been followed even in the simplified context of LQC where, 
however, all the steps could presumably be carried out sufficiently explicitly, see P8_^ for first steps 
into that direction. In full LQG these steps cannot be carried out exactly because the theory is 
too complicated. Thus, in order to make progress one must use suitable approximation schemes. A 
possible starting point could be the following: 

1'. Using kinematical coherent states ipA,E peaked on the constraint surface of the phase space we 
arrive at a subspace of the kinematical Hilbert space TiKin with the property that the constraints 
are satisfied approximately in the sense of expectation values, that is < Cip >~ 0. Denote 
that subspace by 'H'pj^y^ In systems where one can compare predictions by the exact Hphys 
with those by 'H'pf^yg it is often the case that the qualitative answers do not differ much. See e.g. 
|46j . Since there are as many semiclassical kinematical states as classical initial data (classical 
spacetimes), global singularities cannot be determined within this approximate scheme which 
uses the kinematical Hilbert space. 

2'. The formulas for -Pt(t) and Ht are infinite power expansions in the "parameters" t — T |3Uj . 
Hence, close to the singularity, i.e. when r ^ T = (this has to be understood ultimatively 
in terms of operators; we assume that the T, r are chosen such that T = r = corresponds 
to the singularity) we may approximate those formulas by the first few terms and get approx- 
imate expressions P'^(r), H!p that one can hope to handle. These approximations should be 
appropriate in order to answer qualitatively questions about the very early universe. 

^^Alternatively, find an approximation method for the Master Constraint Programme (also using kinematical co- 
herent states and path integral methods, see \17\). 
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3'. The coherent states of item 1' should ideally be chosen as to be preserved by the approximate 
Hamiltonian determined in step 2'. If that is too difficult, peak the kinematical coherent 
states of step 1' on classically singular initial data and then quantum evolve them with the 
(approximate) Hamiltonian Hip, that is r i— >■ exp{iTHp)iljA,E- 

4'. This step is identical to step 4 above with the obvious replacements. 

Notice that one could consider an even easier scheme: Just take any set of kinematical coherent 
states not generated by any (approximate) physical Hamiltonian and consider the unphysical time 
evolution t i— >• ipAo,Eo{'t) '■= '^A{t),E{t) where A{t),E{t) is a one parameter classical gauge evolution 
of the initial data under the flow of the corresponding Hamiltonian constraint (not of any physical 
Hamiltonian)^^ with initial condition A{tQ) = Aq, E{to) = Eq and to ^ tp. Then take any quantity 
P which is not necessarily a (approximate) Dirac observable and study < 4'A{t),E{t)y P4'A{t),E{t) >Kin 
in the limit t ^ corresponding to the classical singularity (in our case A{t) diverges while E{t) 
vanishes at t = 0). This is essentially what we did in the previous section where P played the role 
of the inverse scale factor and the coherent states we used are those of Wile in this case we 
obtained a boundedness result, one may not draw any conclusions from that as the following example 
reveals: 

Consider the Hamiltonian of the hydrogen atom. Let ip^^j^ be the coherent states of the harmonic os- 
cillator peaked at the point (g, p) in phase space. Let g(t),p(t) be the trajectory of a classical electron 
falling radially into the proton from a distance large compared to the atomic radius where t = is the 
point of time when the electron reaches the center. Then it is easy to see that < 4'q{t),p{t), Hip^(^t),p{t) > 
diverges as t — 0. This is because the electron becomes infinitely fast at the center and the expec- 
tation value of the kinetic term diverges (while the expectation value of the potential energy stays 
finite). In contrast, the expectation value < ip^^^it), Hipfj^^{t) >=< ip^^^, Hip^^j^ > is finite where 
'4'q,p{'t) '■= exp(ztff)'?/^g-^p- is the quantum evolved state^^. Thus, while i'q(t),p(t) is nicely peaked on the 
classical trajectory, it does not correspond to the true quantum mechanical time evolution. On the 
other hand, the true quantum mechanical time evolution ipci,p{t) usually does not stay peaked on the 
classical trajectory for sufficiently long time in order to qualify as a semiclassical state representing 
the classical evolution far away from the deep quantum regime closely enough, unless the coherent 
state is actually adapted to the Hamiltonian H. 

Summarizing, ideally the coherent states must be chosen according to the Hamiltonian in ques- 
tion in order to arrive at reliable predictions. If that is too hard, then at least the kinematical 
coherent state should be evolved quantum mechanically but only as long as it remains close to the 
classical trajectory for scales above which the classical theory is reliable. In contrast, the classically 
evolved kinematical semiclassical state not constructed according to the Hamiltonian is inappropriate 
in general which is why the calculation of section E] is promising but inconclusive. 



6 Conclusions 

Concluding, the issue of the presence or absence of the initial singularity in full LQG remains open. 
The results of LQC are important, they are very promising, provide important consistency checks for 
LQG and produce new insights and ideas about LQG. However, as we hope to have demonstrated, 

^^This corresponds to a specific choice of lapse and shift. 
Remarkably, the states ipq^^it) and 4'q{t).p{t) precisely coincide if and only if the underlying dynamics is that of a 
finite or infinite sum of uncoupled harmonic oscillators. 
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they are so far inconclusive for LQG. Therefore sentences such as "LQC is the cosmological sector 
of LQG" or "in LQG the Big Bang is absent" [IHl are misleading for the non expert and should be 
used with due care. 

As our detailed analysis reveals, one must pay attention to the QFT aspects of the full theory 
and cannot rely on the analysis of a quantum mechanical toy model. This is, unfortunately, different 
from the case of the hydrogen atom where experiments confirmed from the outset that the quantum 
mechanical toy model was close to the truth and that the more accurate QED should only give 
corrections. It may still be that LQC is qualitatively right but if that is the case then the mechanism 
for singularity resolution is very different in full LQG and much more elaborate, if LQG is a valid 
candidate for quantum gravity at all. 

In order to resolve the issue about the initial singularity and other cosmological questions in full 
LQG we have proposed an ideal and approximate programme, highlighting the importance of having 
control over the physical Hilbert space and the physical observables. In a zeroth order step of that 
programme we found hints pointing at a resolution of the local initial singularity in LQG by a very 
different mechanism. But it is necessary to carry out all the steps before definite conclusions can be 
reached. Since this programme was not even completed in LQC, again the LQC model could be a 
good starting point to see where the ideal path leads us in a solvable situation. As far as the full 
theory is concerned we will only be able to follow the approximate path and we must learn how to 
control the corresponding approximation errors. 

In summary, the investigation of the comological or black hole sector of LQG, about which the 
LQG inspired LQC model has taught us important first lessons, is an ambitious and fascinating 
research field which will require the joint effort of both cosmologists and quantum geometers in the 
close future. 
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